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Abstract 

Small Dirac masses for neutrinos are natural in models with sin- 
glet fermions in large extra dimensions with quantum gravity scale 
M* ~ 1 — 100 TeV. We study two modifications of the minimal model 
in order to obtain the mass scale relevant for atmospheric neutrino os- 
cillations with at most 0(1) higher-dimensional Yukawa couplings and 
with M* ~ a few TeV. 1) In models with singlet fermions in smaller 
number of extra dimensions than gravity, we find that the effects on 
BR(/i — > &y) and on charged-current universality in ir~ — ► ev, [iv de- 
cays are suppressed as compared to that in the minimal model with 
neutrino and gravity in the same space. 2) If small Dirac masses for 
the singlets are added along with lepton number violating couplings, 
then the mass scales and mixing angles for neutrino oscillations can 
be different from those relevant for u — > e^f and ir~ — > eu, lli>. Thus, 
in both modified models the constraints on M* from BR(^ — ► erf) and 
7r~ — ► eu, [xv decays can be significantly relaxed. Furthermore, con- 
straints from supernova 1987a strongly disfavor oscillations of active 
neutrinos to sterile neutrinos in both the minimal and the modified 
models. 
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1) Introduction 

In models with large extra dimensions and low (~ TeV) quantum gravity 
scale [Q, it is well known that neutrino masses can be naturally small if the 
singlet (right-handed) neutrino propagates in extra dimensions |2|, [| f|. In 
the simplest version of these models, if the quantum gravity scale is smaller 
than ~10 TeV (as motivated by the hierarchy problem) and if the higher- 
dimensional Yukawa couplings are at most 0(1), then it is difficult to obtain 
the neutrino mass scale required for a solution to the atmospheric neutrino 
anomaly. On the other hand, charged- current universality in 7r~ — > eu, [iv 
decays and the limit on BR(/i — > ej) constrain the quantum gravity scale 
to be larger than 10 — 100 TeV, almost independent of the number of extra 
dimensions (assuming large v e — mixing). 

In this paper, we discuss two modifications of this model where it is pos- 
sible to obtain the mass scale required for atmospheric neutrino oscillations 
even with a quantum gravity scale of a few TeV and Yukawa couplings of 
0(1). Furthermore, we will show that in these models, this quantum gravity 
scale can be consistent with charged-current universality in 7r~ — > eu, fii? 
decays and the limit on BR(/i — > e^y). 

Consider for simplicity standard model (SM) gauge singlets Nj = xi) 
in 5.D, where ipi, xi are 2-component Weyl spinors. Assume the following 
terms in the 5D action: 

'-'free = J ^ x< hi iNjT a d a Nj + J d 4 xdy (^fXjNjNj + h.c.) + J d A x iv^o^d^v^ 
5 int = ld 4 x^^^ I (y = 0)h + ^u fXl (y = 0)h + h.cj, 

(1) 

where y denotes the extra dimension assumed to be compactified on a circle 
of radius R and M* is the 5D "fundamental" Planck scale. A, A c are dimen- 
sionless Yukawa couplings in 5D. We have allowed Dirac mass terms (/i) for 
singlets and lepton number violating Yukawa couplings, A c (assigning lepton 
numbers +1 and —1 to Vf and Nj, respectively). The index /i runs over AD 
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while a runs over 5D and T a are the "gamma" matrices in 5D. f is flavor 
index for SM neutrinos and I denotes the index for the singlets and h is the 
Higgs doublet. We have chosen the basis for Ni in which the mass matrix \i 
is diagonal. 

In the effective 4D theory, ipj and xi appear as towers of Kaluza-Klein 
(KK) states, ipi(x,y) = J2 n 1/V%kR e my ^ R ipf l \x) and similarly for Xi{ x iH) 
giving 

5 mass, int = /^E (vi + J^) ^ + 

XfI M * „ ,lS n h M A ^ M * „ 4- h o 1 0\ 

Here we have used the relation, M* +2 R S ~ M Pl (with M Pi ~ 2.4 x 10 18 GeV), 
where 5 is the number of extra spatial dimensions all of which are assumed 
to be of the same size The KK mass term, in/RipjXi\ comes from 
the 5.D kinetic term. Thus, we see that the 4D neutrino Yukawa couplings 
are suppressed by the volume of the extra dimensions or, in other words, by 
the ratio M*/M P i. When h acquires a vev, we get Dirac mass terms for Uf 
with ipj and Xi^\ denoted by m/j ~ XfjvM^/Mpi and rrifj ~ \fjvM*/Mpi, 
respectively (v ~ 250 GeV is the Higgs vev). This analysis goes through for 
any number of extra dimensions. 

Consider the "minimal" model with /i/ = and Xfj = and, to begin 
with, assume one SM neutrino and one N. If mR <C 1, then, to a good 
approximation, u, tp^ and X \ ^ ( n 0) form Dirac fermions with 
masses m and n/R respectively, with mixing between x an d v given by 
~ mR/n <C 1. x decouples and is massless. 

For the case of three SM neutrinos Vf, we can introduce 3 singlets, Ni 
so that the neutrino Dirac mass matrix is (up to small corrections from 

4 Of course, 6 = 1 is not realistic since in that case R is too large. 

5 The number of 2-component Weyl spinors in a Dirac spinor increases with the number 
of dimensions. For example, in 6D, it is 4. Thus, in larger number of extra dimensions, it 
might be possible to give all 3 SM neutrinos Dirac masses using only 1 or 2 singlets. 
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v f -x\ n) mixing g |j) 

{m v ) fI ~ m fI ~ X fi v ^ ~ X fi (^) 10 ~ 4 eV - ( 3 ) 

The v — \ mixing can have significant effects on weak decays to Iv 
as follows. The SM neutrino (weak eigenstate) is dominantly the lightest 
neutrino (with mass m„) with small mixture of heavier neutrinos (with mass 
~n/R): 

(4) 




where v^ 71 ' are the mass eigenstates. The "normalization factor" is 

N ss1 + ?S Lss1 + m!1# =s1 + W (5) 

where we have truncated the KK sum at n/R ~ M*, neglecting an 0(1) 
factor in the summation. For 5 = 2, the sum is log-divergent and we get 
an additional factor of ~ In (Mp^/M*) in the sum. Thus, the decay width 
to z/ - 1 is modified compared to SM since N 2 ^ 1. Whereas decays to z/^'s 
(n ytz o) (if kinematically allowed) are suppressed by small mixing (~ mR/n) 
and have a different phase space. 



For example, consider the decays 7r~ — > ev, /iz/ JTT|. In the SM, T (tt~ — > ev) 
is suppressed by m 2 due to chirality flip. In the extra dimensional scenario, 
the chirality flip can occur on the neutrino instead. In other words, tt~ de- 
cays into ei and v^ ~ ip( n \ i.e., the heavier KK states, through the mvip^ 
term. The large number of KK states up to — m e enhance the effect, 
whereas the same effect (relative to the SM) in the case of n~ — > Hl^^ is 
smaller. For 5 = 2, this effect on T — > ev) gives a lower limit on M* of 
0(1000) TeV for A ~ 0(1) 0. Of course, A ~ 0(1) with M* ~ 1000 TeV 
gives m„ e ~ 0.1 eV, which might be too large. The 7r~ — > eitp^ 1 ' decay width 
scales as |A| 2 (m^/M^) 6 so that for 5 > 3, this effect is smaller than the effect 
of iV 2 on 7r~ decays into e#, and v^ ~ za This modifies the ratio of decay 
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widths to e, fj, since the m's and hence the normalization factors are different 
for v e and the lower limit on M* is O(10) TeV again for — A 2 | ~ 0(1) 
ll|]. We do get the mass scale ~ 10~ 5 (eV) 2 required for a solution to the 



solar neutrino anomaly via matter-induced oscillations [15|] for M* ~ 0(10) 

TeV and A ~ 0(1). 

Coherent conversion of SM neutrinos to singlet neutrinos (due to the 

above mixing) in a supernova (SN) results in energy loss, reducing its active 

neutrino flux. Since the SM neutrinos (unlike the singlet neutrino) have weak 

interactions with the matter in the SN core, these oscillations are enhanced 

by the MSW effect. These resonant oscillations are possible only if the mass 

of the sterile neutrino state is not larger than ~ V EV ~ 10 keV, where 

E ~ 100 MeV is the neutrino energy in a SN and V ~ 10 eV is the potential 

in its core |J. The survival probability of SM neutrino can be approximated 

by the product of survival probabilities in each resonance "crossed" by the 

SM neutrino as it travels out of the SN core, P vv ~ Yl n P n H Hi- P n is 

th 

approximately independent of the mass of the n resonance and is given by 
1, 1 



/ 7r4mV CO re\ / mf \ la , 

Pn W eXP {- 2 — B— J ~ GXP {-WHeW J ' (6) 
where r CO re ~ 10 km is the radius of the SN core. To explain the atmospheric 
neutrino anomaly via oscillations, we require m? (for z/ M or u T ) ~ Am^ m ~ 
10 -3 (eV) 2 [^] so that each P n is ~ 1/3. Thus, even if only one resonance is 
crossed, the energy loss from a SN due to sterile neutrinos will be comparable 
to that due to active neutrinos. The agreement of the measured duration and 
number of neutrino events from SN1987a with the prediction of SN models 
(taking into account theoretical uncertainties) indicates that a new channel 
for energy loss from the SN should be less effective than the standard neutrino 
channel; otherwise the neutrino signal duration will be halved, which is not 
allowed |TjJ. Therefore, in this case, the measurement of SN1987a neutrino 
flux implies that no resonance can be crossed so that the mass of the lightest 
KK state, which is 1/R, should be larger than ~ 10 keV ||. This implies 
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that the KK states are much heavier than SM neutrinos and hence both the 
solar and atmospheric neutrino anomalies have to be explained by oscillations 
among active neutrinos, governed by the mass matrix m v above. 

Loop diagrams involving KK neutrino tower and longitudinal W con- 



tribute to — > e7 [0, §] and (g — 2)^ |9|, [10(1 . The coefficient of the dimension-5 
operator relevant for these two processes, F^la^J (I = e, fi, r), generated by 
the KK neutrino exchange, is approximately 

mi j- ^ R 2 mi mm^ M\ x 

The amplitude is enhanced as compared to 4.D case by the large number of 
KK states (we have truncated the KK sum at M*). Since the neutrino masses 
are given by m (see Eq. (|3|)), it is obvious that (as in AD) there is a direct 
correlation between neutrino oscillations and the contributions to /j, — > e7 
from these loop diagrams. For example, in the two flavor case, we get 



(mm^) 
v / < 



M 
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where m Ul , m U2 and # are the masses and mixing angle of the two neutrinos 
obtained from Eq. (|^). This was used in ]7| to obtain lower limits on M* 
of O(10 - 100) TeV for 9 ~ vr/4 and \m\ - m 2 U2 \ ~ 10~ 5 (eV) 2 (as relevant 
for large mixing angle solar oscillations). In the case of three flavors with 
2 — 3 mixing angle 0, but no 1 — 3 mixing, the above expression is simply 
multiplied by cos0. 

In the minimal model, we see from Eq. (^) that to get m 2 , ~ Am^ m we 
require M* ~ 100 TeV if A ~ 0(1). Such high values of M* are disfavored by 
the motivation to solve the hierarchy problem [0J. Of course, we can choose 
A > 1 and obtain m 2 u ~ Am 2 ltm for M*~10 TeV Q, but then the (A + 5)D 

6 For Af» ~ 10 TeV, it is not possible to obtain Ar7i| tm even for A ^> 1 since, due to 
the normalization factor, there is an upper limit on m u , for a given Af» . 
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theory might reach strong coupling. In other words, from the AD point of 
view, the Yukawa coupling ~ XM^/Mpi ~ m u /v is very small at tree level (or 
at low energies) even though A ^> 1 in this case. But, since the 4D coupling 
"runs" with power of energy due to the multiplicity of KK states, it might 
reach its Landau pole near M*. For this reason, we will consider A ~ 0(1) 
throughout this paper. 

With the motivation of obtaining the neutrino mass scale Am^ m with 
M* ~ TeV and A ~ 0(1), we now study two modifications of the minimal 
model. In the first model, the singlet neutrino propagates in a sub-space 
of the full extra dimensional space where gravity propagates. In the second 
model, we consider the effect of non-zero Dirac masses for the singlets and 
of lepton number violating couplings, A c . We will also keep an eye on the 
correlation between neutrino masses and contribution to fi — > ej and the 
effect on 7r~ — > eu, fiv in these models. 

2) Sub-space 

First, consider the case where singlet neutrino propagates in 5 U < S di- 
mensions @. Assuming that all extra dimensions are of size R, in this case, 
we get the neutrino Dirac mass matrix 

A// . /M* \ s » /s 

V ~ AfiV [ — 

{RM*) 5u KM pi- 

Thus, for 5 U = 5 and 5 = 6 and with M* ~ TeV, A ~ 0(1), we get m 2 ~ 
Am^ m 0. To obtain the neutrino mass scale required for a solution to 
the solar neutrino anomaly via matter-induced oscillations, Am^j ~ 10~ 5 
(eV) 2 , we can choose the corresponding A ~ 0(0. 1). 

In this case, the dimension-5 operator relevant for fi — > ej has the coeffi- 
cient 



(m v ) fI ~m fI ~ ■ = v~\ fI v(— - . (9) 



n 2 v 2 M 2 ^ n 2 g n 2 v 2 M 2 



I 2 2 

m:. — mf 





Su/5 

1 






(Ml, 
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where vn Vx , m U2 and 6 are the neutrino masses and mixing angle as obtained 
from Eq. @ . We see that this contribution is suppressed compared to 5 U = 5 
(for the same values of m and M*) due to smaller number of KK states. Thus, 
the lower limits on M* of O(10 — 100) TeV obtained for the minimal model 
[0] (for 6 ~ 7r/4 and |m 2 — mL\ ~ 10~ 5 (eV) 2 ) can be relaxed by a factor of 
O(20), assuming 5 V = 5 and 5 = 6. The lower limit now becomes M* ~ few 
TeV. 

In terms of A (instead of m), the coefficient of the dimension-5 operator 

is 

™,( AAt L , , 

8tt 2 M, 2 ' 1 J 

i.e., for fixed A, it is independent of 5 U or 5 |§. But, with 5 U = 5, 5 = 6 
and M* ~ few TeV, we require ( AAH ~ O(10 -2 ) to obtain the mass scale 
for solar neutrino oscillations (as shown above). Whereas with 5 U = 5 and 
for M* ~ O(10 - 100) TeV, we require fAA tN ) ~ 0(1). Hence the above 
coefficient is the same for these two parameter sets (which give the same 
neutrino masses), in agreement with the analysis in terms of m. 

Reference || also considers the constraints from \x — > 3e and /j, — > e con- 
version in nuclei. For these processes, loop contribution due to KK neutrino 
tower to the effective Z-fi-e coupling (in addition to the 7-/x-e coupling) has 
to be included. This coupling depends on ~ A 4 t> 2 /M 2 || (dropping the fla- 
vor indices on A for simplicity), unlike the 7-ju-e coupling which depends on 
~ A 2 f 2 /M 2 as above. The experimental bounds on \i — > 3e and fi — > e con- 
version in nuclei give the lower limit M*/A 2 ~ 200 - 300 TeV [| which, for 
A ~ 1, is stronger than that from \i — > erf. However, in this sub-space case, 
as mentioned above, to get m 2 ~ Am 2 Q j with M* ~ few TeV, we require 
A ~ O(0.1) and hence this quantum gravity scale is consistent with /x — > 3e 
and fi — > e conversion in nuclei. 

In this sub-space scenario, the effect on 7r~ — > ev, \w decays is due 
to the normalization factor A^ 2 (since 5 V = 5) and thus also depends on 
J2 n { m R) 2 /n 2 ( see Eq. (|5])). As above, this factor is smaller than in the 
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minimal model for the same values of to and M*. Therefore, the lower limit 
on M* from 7r~ — > eu, \w decays is also reduced from O(10) TeV (obtained 
for the minimal model ||11|| , assuming A ~ 0(1) which gives m 2 ~ 10~ 5 (eV) 2 ) 
to ~ few TeV. 

With 5 = 6, we get 1/R ~ 0(10 - 100) MeV so that the SN1987a 
constraint is satisfied. 
3) See-saw 

Next, consider the case with non-zero Dirac mass /i/(<C M*) for the sin- 



glets p, [12|, [L3| . In the limit /i 3> m and m c »s 0, the linear combina- 



tions ~ fxj^ + rrifi/ fij + n 2 / R 2 u^j form Dirac pairs with ^j^'s of masses 



' /if + n 2 / i? 2 . Of course, with A c ~ 0, there are 3 massless neutrinos which 
are dominantly f/s, even though a priori all fermions have mass terms. The 
reason is that in this case we can define a conserved lepton number with 
charges +1 for v and x an d —1 for ip so that only Dirac masses are allowed. 
Thus, there are 3 "unpaired" fermions with charge +1 which are the z/j's. In 
the minimal model (fi = 0), these massless fermions are xf^ s as mentioned 
earlier. We require lepton number violation, for example, m c ^ 0, so that 
these massless neutrinos can get Majorana masses VfVp |I2, [13]| . Then, the 
see-saw mechanism (see Fig. P gives 



to,. 



7/' 



i 



/xj + n 2 /R 2 



In this case, one singlet N suffices to give masses to all 3 SM neutrinos, but 
in general, one can have many singlets with different /i/'s. 

In the case of 5 > 2, the see-saw is "divergent" due to J2 n ^/ n2 - If we 
truncate the KK sum at M*, then we get the neutrino Majorana mass matrix 



'//' 



(mam c T + rrfam , 
(A,A- + A>A-) ff ,^. (13) 
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Figure 1: See-saw mechanism for generating SM neutrino masses in the case 
Hi 7^ and m c ^ 0. There is another diagram obtained by m «-» m c , ?/>j n ' ) <-> 

(n) 

X/ • 

The KK sum is log-divergent for 5 = 2 so that the above result is multiplied 
by a factor of ~ In M 2 / (jj 2 + l/R 2 ) ~ O(10) and was mentioned in f|TJ]. 
We see from Eq. (fL^D that for M* ~ TeV, <5 > 2 and A, A c ~ 0(1), we get 

~ Am^ m by choosing /i ~ eV. If all singlets have Dirac masses of O(eV), 
then to get m 2 u ~ Am^j, we can choose the corresponding A, A c ~ O(O.l). 

We next consider the constraints from SN1987a on this model. The anal- 
ysis is different than in the minimal model since in this case m u ^ m (the 
latter governs the mixing between SM and sterile neutrinos and hence the 
energy loss in SN1987a) and also the masses of the sterile neutrinos are given 



by ~ J jj 2 + n 2 /R 2 instead of n/R. Thus, one way to evade the SN1987a 
constraint, independent of the number and the size of the extra dimension, 
is to choose fii ~ 10 keV so that no resonance can be crossed in the SN1987a 
core. Then, the solar and atmospheric neutrino oscillations have to involve 
only active neutrinos. With fi ~ 10 keV, we require AA C ~ O(10 -4 ), O(10~ 5 ) 
to get ml ~ Amj^. m , Am^j, respectively, i.e, A, A c ~ 0(1) will give too 
large m u (see Eq. ([13])). 

If 8 > 4, then we get 1/R> 10 keV so that again no resonance is crossed 
and the SN1987a constraint is satisfied for any value of fi. 

The remaining cases are fii 10 keV and 5 = 2, 3. The number of 
resonances crossed, i.e., the number of sterile neutrino states lighter than 
10 keV, is n r es ~ (R 10 keV)* 5 , and we have to consider each value of S 
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separately. 

For 5 = 2 and M* ~ 1-10 TeV, we get l/R ~ 0.01 eV so that n res ~ 10 12 . 
As mentioned earlier, we require 

Puu « (Pn)" ros ~ 1/2 so that the SN energy 
loss due to sterile neutrinos is less than that due to active neutrinos. Thus, we 
get the constraint m 2 , m c2 ~ 10~ 15 (eV) 2 . Then, the active neutrino masses 
are too small to account for solar and atmospheric neutrino oscillations. The 
v e — ^( n )/-?/>( n ) mixing is also too small to be relevant for explaining the solar 
neutrino anomaly by oscillations of v e to sterile neutrinos. Thus, /x <C 10 
keV is ruled out by SN1987a constraint. 

For S = 3 and M* ~ 1 - 10 TeV, we get l/R ~ O(100 eV - 1 keV) so that 
77, res ~ 10 3 — 10 6 0. Hence, the SN1987a constraint, P uv ~ 1, requires m 2 , 
m c2 ~ 10 -6 — 10" 9 (eV) 2 , i.e., A, A c ~ 0(1) might be allowed (depending on 
M*). Thus, we can choose /x ~ eV to get m 2 ~ Am^ m for M* ~ TeV and 
A, A c ~ 0(1) while (marginally) satisfying the SN1987a constraint. 

Of course, for 5 > 3 and for any value of /x, the sterile neutrinos are heavier 
than l/R ~ 100 eV so that they cannot be directly involved in oscillations 
of SM neutrinos in the sun or in the atmosphere. 

The coefficient of the dimension-5 operator, F^la^l, generated by ex- 
change of KK neutrinos, is similar to the earlier case, with an additional 
contribution from m c : 

mi mm* + m c m c * Mp t 

Thus, we see that both (g — 2) M and BR(/i — > erj) depend on (mm* + m c m c *^j 
whereas m v depends on (m /x m cT + m c /x m T ^j. Therefore, it is clear that 
the parameters for neutrino oscillations and for /x — > &y, (g — 2)^ may not be 
related in this class of models. In particular 

A (/x — > ej) oc (mm^ + m c m c ^ ^ \ m u! ~~ m u 2 1 s ^ n 0-**) 
"' R and hence nres is very sensitive to the value of M* . 
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where the masses m vi , m V2 and mixing angle 9 (relevant for v e -v^ oscillations) 
are obtained from m u in Eq. (|13|) . 

To illustrate the above point, consider for simplicity only v e and and 
suppose we have 2 singlets with the same /i's. Assume 

(16) 
(17) 

with p, a ~ 0(1). In this example, p — > ey depends on 

(18) 

\u it o~ -t- p~ i 
whereas 

m v oc f ° n P ] (19) 





so that z/ /t — z/ e mixing can be large (with non-degenerate neutrinos) as re- 
quired for solar oscillations with large mixing angle, but the loop contribution 
to p — > ey is zero. 

Even if the flavor structures of m and m c are similar, BR(/i — > cy) might 
not be correlated with neutrino masses since, for given M*, p — > ey depends 
only on A and A c , whereas m u depends also on p. For 5 = 2, as explained 
above, to satisfy the constraints from SN1987a, we might have to choose 
p ~ 10 keV and hence AA C ~ 10~ 5 to get m 2 u ~ Am^j. Then, ~ few TeV 
is consistent with the limit on BR(/i — > ey). For 5 > 3, p is not constrained 
by SN1987a. So, we can choose the relevant A, A c small enough such that 
M* ~ few TeV is consistent with p — > ey and, at the same time, we can get 
the required m Ve , m u ^ by choosing the corresponding p^s appropriately. As 
mentioned in section 2, loop contributions to p —>■ 3e and p — > e conversion 
in nuclei also depend on A and A c and hence are suppressed if A, A c are small. 
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A similar analysis shows that M* ~ few TeV can be consistent with 
7r~ — > eu, \xv since, for given M* and for /i <C m w , the effect on tt~ decays 
also depends only on A and A c . 

4) Conclusion 

In summary, we have studied two "non-minimal" models with singlet 
neutrino in large extra dimensions and TeV scale quantum gravity. These 
models can accommodate the mass scale relevant for atmospheric neutrino 
oscillations even with 0(1) higher-dimensional Yukawa couplings and ~ 
a few TeV, unlike the minimal model. The first model has singlet neutrino 
propagating in smaller number of extra dimensions as compared to gravity 
whereas the second model has small Dirac mass terms for the singlets and 
lepton number violating couplings. In both models, the constraints on M* 
from BR(/x — > c-f) and tt — > eu, fiu decays can be significantly weakened as 
compared to the minimal model so that ~ few TeV is consistent with 
these decays. Also, due to the SN1987a constraint, active-sterile neutrino 
oscillations are strongly disfavored in both the minimal and the modified 
models. 
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